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Abstract Using the expression of the dynamical gluon mass obtained through 
the operator product expansion we discuss the relevance of gluon mass effects in 
the decays V — > hadrons (V = J/ip, T). Relativistic and radiative corrections are 
also introduced to calculate a s (m c ) and a s (mb) comparing them with other values 
available in the literature. The effects of dynamical gluon masses are negligible for 
T decay but important for J/ip decay. 
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Since the original suggestion by Appelquist and Politzer U, the annihilation of 
heavy quark pairs into gluons or gluons plus a photon has been recognized as an 
apparently excellent process for testing the basic ideas of QCD and for measuring 
the strong coupling a s with high precision. 

One can determine a s at low energies by comparing the experimental T and J/i/j 
branching ratios 

t~> r(V -> ggg) 

Rv = T(V -> ee) ' V = T > J ^ C 1 ) 
with the theoretical prediction. This is given, at the lowest order in perturbation 
theory, by 

R 10 (?r2 ~ 9) n 3 (2) 
Rv ~ Slnel a(M v ) ' (2) 

where e g is the quark charge, a is the electromagnetic coupling and My is the 
quarkonium mass. In addition, the experimental data for the radiative process 
V — > 7 + X is an important source of information about the QCD dynamics, which 
predicts a photon spectrum nearly linear in z = 2E 1 /My, where E 1 is the photon 
energy. 

The match of the experimental and theoretical determinations of a s and the 
photon spectrum with high accuracy has originated a long series of discussions since 
the pioneer work of Ref . • Corrections due to relativistic and radiative QCD ef- 
fects were shown to modify Eq.(^). Other effects, as the existence of an effective 
gluon mass or beyond standard model contributions, are also possible explanations 
for the claimed differences between theory and experimental data. In this work we 
introduce the concept of a dynamical gluon mass and verify how it affects quarko- 
nia decays. The fact that we consider a momentum-dependent gluon mass clearly 
modifies previous discussions on this problem. 

Parisi and Petronzio were the first to discuss some discrepancies between 
QCD results and the ones obtained from J/ifj decays. They introduced an effective 
gluon mass to improve the comparison between experiment and theory in the case 
of J/ip decays. This is not the unique alternative used to solve the problems in this 
context but this line of thought does give reasonable results, and the possible effects 
of a gluon mass in quarkonia decays were considered several times in the literature 
(see |§-|| and references therein). In Refs. 0, |3], |J the value of the gluon mass 
was determined calculating the CM energy spectrum of direct photons produced by 
J/ip and T decays (J/ip,T — > ggj). Afterwards the spectra were compared with 
experimental data, which clearly favoured the case of a massive gluon and fixed its 
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mass. It is interesting to note that the explanation of the photon spectrum in the 
J/ip decay demanded a gluon mass M g ~ 0.66 GeV, whereas in the T decay the 
most suitable mass value raised to M g ~ 1.2 GeV 0, |]. There are some points 
that may be criticized in these results. With the knowledge accumulated in the past 
decades about the theory of the strong interactions it is obvious that we cannot have 
a bare gluon mass in QCD. In Refs. [f] the gluon masses differ in the J/tp and 
T decays in an unnatural way: as we increase the quarkonium mass the gluon mass 
needed to explain the photon spectrum is also increased. Finally, the explanation 
of the bending of the photon spectrum at large z in these decays may be due to the 
radiation of additional gluons as described by Field 0. 

It is clear that QCD does not admit a bare gluon mass. However, this the- 
ory may have a dynamical gluon mass M and recent simulations of QCD on the 
lattice strongly support this possibility ||. Therefore, the proposal of Parisi and 
Petronzio may indeed be realized in Nature, but in a more subtle way. This will 
be our working hypothesis. The next question is how to consider this dynamical 
mass in quarkonium decays. The introduction of a gluon mass scale in the calcula- 
tion of some hadronic processes has been performed in a rather heuristic way (see, 
for instance, Ref. [fnj). Only recently a more formal handling of gluon masses has 
been discussed [JTTJ]. In principle, we can justify the approach of Ref. [[H]] within the 



dynamical perturbation theory proposed by Pagels and Stokar many years ago [12 
which basically imply in the use of the running gluon mass in the gluon propagator. 
This is the approach that we will follow here. 

In order to obtain the asymptotic behavior of the gluon mass we can rely on 
the operator product expansion (OPE). This asymptotic behavior was obtained in 



Ref. 13 



where P 2 (= —p 2 ) is the gluon momentum in Euclidean space, {^■G til 'G is the 
gluon condensate fllH , and iV = 3 for QCD. At this point we see the difference with 
previous work || f|]: the dynamical mass is connected with the gluon condensate 
and it decreases with energy. Therefore, we can expect different results from the 
ones of Refs. and 0. 

Equation (||) gives the asymptotic behavior of the gluon mass but we also need 
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its expression in the infrared region. If we define 



("9(iV 2 



1/2 



u , (4) 



it was recently verified that the value of the dynamical gluon mass in the infrared 
is well described by the OPE value frozen at the scale m g [15|, i.e. we can write 

4 

777 

M 2 g (P 2 ) = mpiml - P 2 ) + -pf 9(P* - m 2 ). (5) 

This is the expression which will be used in the calculation of quarkonium decays. 

As discussed in Ref . M , the gluon mass implies that the branching ratio Ry must 
be changed to Ry ■ fziji), where 

, / x T(V ^ ggg)\ mg fR . 
r(V -> ggg)\ mg =o 

is a function of 77 = 2m g /Mv- This simple factorization is a consequence of the fact 
that most of the gluon mass contribution in the calculation of T(V — > ggg) comes 
from the phase space integration |2], [J. Therefore, we proceed in the same way as in 
H to compute ^3, with the difference that the mass is now a function of each final 
state gluon momentum. 

In order to compute the function / 3 as a function of 77 we made use of the 
package COMPHEP | 16|| . Initially we considered the case where a bare gluon mass 
was introduced in the matrix element and in the phase space. The result is given by 
the solid line in Fig.l, which is identical to the one of Ref.p[. We then introduced a 
bare mass only in the phase space integration, obtaining the thin dashed line of Fig.l 
(compatible with Ref.[[|). This result shows that the mass effect can be accounted 
for within this simple approximation. Once we verified that this approximation 
reproduces most of the mass effect, the calculation to obtain _f 3 was performed by 
Monte Carlo integration considering the dynamical gluon mass contribution (Eq. 
(^D) only in the phase space. The result is the thick dashed line depicted in Fig.l. 

As discussed before, in Refs. and 0, the value of the gluon mass was obtained 
by calculating the CM energy spectrum of direct photons produced by J /if) and T 
decays (J /if) , T — > ggj), yielding very different values in each case. With the gluon 
mass obtained from the fit of the radiative quarkonia decays, the value of / 3 was 
determined from a curve similar to Fig.l, leading to a new determination of a s . 
Here the gluon mass m g is given by Eq. (f|). If we consider the gluon condensate 
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Figure 1: Function f^{rj) representing the gluon mass corrections. The solid line 
represents these corrections when the effect of a bare gluon mass is taken into account 
both in the matrix element and phase space. The thin dashed line is obtained 
considering the effect of a bare gluon mass only in the phase space. The thick 
dashed line is our result obtained by taking into account a dynamical mass as given 
by Eq. (|). 



value given in Ref. 0, [r7| , and the gluon mass obtained in Refs. p, [L(J, [L8|, we can 
assume m g « (0.64±0.20)GeV, where the uncertainty is large enough to express the 
various limitations present in the dynamical gluon mass estimates. It is interesting 
to see that the central value is of the order of the value determined by Consoli and 
Field || U to explain the J/ip — > X + 7 decay (m g m 0.66 GeV). With this value 
for the gluon mass we certainly do not explain the direct photon spectrum in the 
inclusive decay of the T. However, the explanation of this process may depend on a 
complex relationship between radiation of additional gluons as described by Field |?J 
and the existence of a dynamical gluon mass. 

As one can notice in Fig.l, the new curve (the thick dashed one) is flatter than 
the curves of Refs. [§, 0] for small values of rj. This is not so surprising since 
as m g decreases f$ tends to 1. But this feature is strengthened in the case of a 
dynamical mass whose momentum dependence is given by Eq.(|5]), because for small 
T) (which means small gluon mass or large quarkonium mass) the integrated phase 
space is larger than in the case of a bare gluon mass. This is so because we have 
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configurations where at least one of the gluons resulting from the quarkonium decay 
has its mass effect sharply cut-off by a large momentum. 

Using the values of f 3 from Fig.l we can now calculate a s {m c ) and a s (m&) . 
Assuming m g = (0.64 ± 0.20) GeV, we have (rj = 2m g /M v ) 

= (0.41 ±0.13), (7) 

and 

r/ T) = (0.14 ±0.04), (8) 

and from the curve of Fig.l, we obtain 73(77), i.e. the factors due to gluon mass 
corrections are 

fafaVM) = 0.47 ± 0.30 , f 3 (r] {T) ) = 0.94 ± 0.03. (9) 

Our value of fz{rf J t' 1 ^) is consistent with the one in Ref.Q, £|, but f^irf^) is ap- 
proximately 30% larger and almost identical to 1. If the mass was also introduced 
in the matrix element, the values of f% would be increased by a few percent and no 
signal of the gluon mass would appear for f^{r]^). 

We can now compute a s from the J/ip and T decays, but to do so we have to 
take also into account the relativistic and QCD corrections. The QCD corrections 



are introduced through the factor (see p"9|l) 



a-.s 



7T 



by) (10) 



where by = 1-6, 0.43 for J/ip and T, respectively. 

Equation (^) is obtained assuming that the qq pair annihilation occurs at a 
point. Relativistic corrections arise if we consider that this process occurs over a 
finite volume of radius ~ l/m q . Several papers have dealt with this issue. We 
follow Ref. [[2(J and use their correction factor: 7^ = (0.31 ± 0.03) for the J/ip and 
T W = (0.69 ± 0.07) for the T. 

Equation (0) rewritten with the correction factors mentioned above becomes 



m,, \ 10 [~- - !)) ,, 



R v = fM 1 + ^ 7 ai 2 (M \ « • (11) 
7r / 8l7re^ a(My) 



We then use the experimental values of Ry from plj and 0] to calculate a s {m c ) 
and a s (mb): Rj/^ = 10.1 ± 0.9 and Ry = 32.6 ± 0.8. Our results are 

a s (m c ) = 0.35 ±0.07, (12) 



6 



and 

a s (m b ) = 0.206 ±0.008. (13) 

The values of Eqs. (|12|) and ([il|) are compatible with the ones of Refs. 0, ^(| . Note 
that Ref. |3J] does not consider the relativistic corrections of Ref. [P0[| . On the 
other hand, our result for ^(t/^) is about 1.6 times larger than the one of Refs.0, 
|J]. We find that the a s (m(,) determination through T decays is barely affected by 
the dynamical gluon mass effects and, considering the relativistic corrections, it is 
totally consistent with the one of Ref. || (where the relativistic corrections were not 
considered). Contrary to the values of a s {m c ) obtained in the references quoted 
above our result has a large error, and most of it, as we discuss next, is due to the 
uncertainty in the gluon mass. 

Assuming that the gluon propagator has an infrared mass scale f|, we see sev- 
eral limitations in order to determine a s with high confidence level through heavy 
quarkonium decays. The main problem of this approach is that the dynamical gluon 
mass scale is poorly known. The several determinations of this mass scale are char- 
acterized by a large range of possible values |8|, [It], I8"f . Even if we consider that this 



mass follows the behavior predicted by OPE |13|, we know that the gluon conden- 



sate in opposition to the fermionic one, is known with a large uncertainty fl7 |. The 
extrapolation of this value to the infrared region also involves a series of approxima- 



tions |15| . Finally, if we attempt to verify the consistency between measurements of 
a s at different mass scales, we also may be in trouble if we do not take into account 



the gluon mass effect in the running coupling constant [[22]|. 

If the gluon propagator has a dynamical mass scale the strong running coupling 
constant is modified, at low energy, to f| |22| 

slV ' " (lliV - 2n f ) ln[(Q 2 + £m 2 g )/A 2 } ' 1 ' 

where nf is the number of fermion flavors, and £ ~ 4 (or even larger pIJ) is a 
parameter determined in Ref. ||. This behavior of the running coupling constant 
with respect to the gluon mass appears when solving the Schwinger-Dyson equations 



for the gluon propagator and the trilinear vertex. The determination of Eq. (|14"|) and 
its higher order corrections is a much more complex problem than the determination 
of the perturbative expression of a s (Q 2 ), and, unfortunately, is not under complete 
control as the perturbative one. Therefore, at the scale rrib and up we do not expect 
any substantial difference in the phenomenology if we use the perturbative running 
coupling constant when checking the consistency of a s (ml), obtained through Eq.(Ef), 
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with others measurements of a s . At the scale m c this is not true anymore. In this 



case Eq.([L4]) and its higher order corrections (not yet computed) have to be used. 
At leading order of both coupling constants (perturbative and nonpertubative), the 
difference at the scale m c in using one or another expression amounts a factor roughly 
given by In {m 2 c + ^m 2 g ) / \n.m 2 c « 1.6. This number is large enough to interfere in 
high precision measurements of the running coupling constant using J/ijj decays, 
and a better understanding of the gluon mass dependence of the running coupling 
constant is needed. 

In conclusion, we have presented an analysis similar to that of Refs. ||, H f|]. 
Notwithstanding, we considered the effects of a dynamical gluon mass, as given 
by Eq. @, which is compatible with the asymptotic behavior predicted by OPE 
[Oj. Monte Carlo integration of the phase space in quarkonium decays provided a 
new curve for the parameter fsfy), which is flatter for small values of rj than the 
previous ones, and reflects the momentum dependence of the dynamical gluon mass. 
Assuming m g « (0.64 ± 0.20)GeV, we have obtained new values for fs(r] J ^) and 
/ 3 (?7 T ). With these values we calculated a s (m c ) and a s (mb) including QCD as well 
as relativistic corrections. The values of a s which we have obtained, within the 
assumed error bars for the gluon masses, are compatible with the ones in Refs. 
2(J. Contrarily to previous analysis we verified that a s {m h ) is not affected by the 



existence of a dynamical gluon mass. More rigorous determinations of the gluon 
mass scale are necessary in order to determine a s (m c ) with high precision by using 
heavy quarkonium decays. 
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